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ABSTRACT 

We discuss in the framework of the excursion set formalism a recent discovery from 
N-body simulations that the clustering of haloes of given mass depends on their forma- 
tion history. We review why the standard implementation of this formalism is unable 
to explain such dependencies, and we show that this can, in principle, be rectified by 
implementing in full an ellipsoidal collapse model where collapse depends not only 
on the overdensity but also on the shape of the initial density field. We also present 
an alternative remedy for this deficiency, namely the inclusion of collapse barriers for 
pancakes and filaments, together with the assumption that formation history depends 
on when these barriers are crossed. We implement both these extensions in a gener- 
alised excursion set method, and run large Monte Carlo realisations to quantify the 
effects. Our results suggest that effects as large as those found in simulations can only 
arise in the excursion set formalism if the formation history of a halo does indeed de- 
pend on the size of its progenitor filaments and pancakes. We also present conditional 
distributions of progenitor pancakes and filaments for low-mass haloes identified at 
present epoch, and discuss a recent claim by Mo et.al. that most low-mass haloes were 
embedded in massive pancakes at z ~ 2. 

Key words: methods: statistical - cosmology: theory - galaxies: clustering - galaxies: 
haloes - large-scale structure of Universe 



1 INTRODUCTION 

The excurs i on se t formal is m l|P ress fc SchechteJ 
[ Bond et alJ Il99lt iBowerl Il99ll : ILacev fc Cold Il993: 
iMo fc Whitd [l996) has been remarkably successful in 
describing several characteristics of the halo population 
in numerical simulations of structure formation, such as 
their unconditional and conditional mass functions, their 
merger rates, the halo bias and more. In general terms 
the approach implies picking a "particle" at random and 
smoothing the linear density field over ever smaller spheres 
around it, until the criterion for collapse at some redshift 
z is satisfied. The mass in the sphere is then identified as 
that of the collapsed object to which the particle belongs. 
This criterion is typically a constraint on the linear density 
contrast obtained from a collapse model, either spherical or 
ellipsoidal. The spherical collapse model gives a simple crit- 
ical value for collapse S c « 1.686 a bove which a patch is sai d 
to have collapsed at present time JPress fc Schechterlll974) . 
For the most popular implementations of the ellipsoidal 
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model this is modified only by making the collapse criterion 
scale-dependent 8 c (a(R)), because shape parameters are 
approximated by their exp ectation values iSheth et alJ 
l200lt ISheth fc Tormenll200a) . Thus determining whether a 
patch has collapsed or not requires knowledge only of the 
density contrast 8 and the smoothing scale R. 

When smoothing the density field a specific choice of 
filter shape must be used. An intuitively appealing choice is 
the top-hat filter for which both visualisation of the prob- 
lem and the assignment of a mass to each collapsed region 
are straightforward. However, since the fc— modes of the lin- 
ear perturbation field are independent, it is mathematically 
more convenient to use a sharp fc— space filter, since varying 
resolution then corresponds to a Markov random walk in the 
density contrast. As it turns out, this is also the filter which 
best reprod uces the shape of mass functions seen in N-body 
simulations JPercivaJl2001^ . although a different filter then 
needs to be employed to assign a halo mass to each smooth- 
ing scale. It is surprising that this simple formalism, based 
only on the linear density field, is so successful in explaining 
many of the features seen in N-body simulations. 

Despite its successes the approach clearly has shortcom- 
ings. Perhaps the most studied arises from the inapplicabil- 
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ity in detail of the physical model which underlies it. This 
stems from the fact that few particles actually lie at peaks 
of the smoothed initial density field. Most lie somewhere on 
the peripheral slopes of the peak associated with the halo 
they will actually belong to at the target redshift. As a re- 
sult, there is a relatively poor match between the halo mass 
predicted for individual simulation particles by the random 
walk associated with the linear density field and the mass 
of the actual halos in which they find themselves as a result 
of the nonlinear evolution llBond et alJll99lt IWhiteill996t 
ISheth et alJl200lTl . The excellent agreement of the statis- 
tics of simulated halos with excursion set theory predictions 
thus does not e xtend to good object-by-object agreement. 
iBond fc Mversl il996l) devised an algorithm for explicitly 
taking account of the non-locality of the collapse process 
by assigning a surrounding "patch" to each "peak" in the 
initial conditions. This adds considerable complexity, how- 
ever, and as a result is rarely used. Despite this problem, 
the excursion set approach has proved remarkably useful. 

Another weakness which has received recent at tention 
JWang et alJl200fJ; ISheth fc TormerJ 12004 see also IWhitd 
l^9w is the fact that in its simplest and most used approx- 
imation, the approach does not allow correlations between 
halo properties defined by the character of the random walks 
on opposite sides of the barrier crossing. One example where 
such corre l ations are seen in simulations is the discovery by 
iGao et alJ l l200Sh (GSW05) that old haloes are more clus- 
tered than young haloes of similar mass, or, more generally, 
that the clustering of halos d epends on their form ation his- 
tory as well as on their mass llCroton et al . 2006, CGW06). 
This is also manifest as a dependence of clustering on prop- 
erties such as halo c oncentration or substructure fraction 
iWechsler et al . 2006). All these effects require some degree 
of correlation between environment density and halo for- 
mation history. In the standard 1-dimensional excursion set 
approach such correlations are impossible by construction. 

There are two ways to extend the excursion set ap- 
proach to address these effects. For a m ultidimensional 
random walk of the kind pioneered by IChiueh fc Led 
(2001)(CL01) the extra variables can carry information 
across the barrier. A particle with a high environment den- 
sity can cross at a different point in terms of these extra 
variables than a low density particle, and this can then be 
reflected in the random walk at smaller scales, hence in the 
formation history of the particle's halo. A second way to ac- 
count for such correlations arises if the "shape" of a halo's 
formation history depends on different barriers than does its 
mass. For example, if we allow for the possibility that the 
assembly history of a halo could be significantly impacted by 
the size of its progenitor pancakes and filaments, an interest- 
ing scenario appears: it is then possible that the progenitor 
pancake scale is larger than the present halo scale, meaning 
the pancake barrier crossing lies between the halo collapse 
crossing and the scale defining the environment. This makes 
possible some correlation of formation history with environ- 
ment density for halos of given mass, as seen in simulations. 
In this paper we will study both these possibilities in detail. 

In order to allow for multi-scale correlations we imple- 
ment the 6D Markov approach pioneered by CL01 using re- 
alisations of the deformation tensor rather than merely the 
density contrast, thereby utilizing the full shape informa- 
tion. We introduce a new set of shape parameters which al- 



low us to construct a halo collapse barrier w hich reproduces 
very w ell the dynamical collapse model of IBond fc Mversl 
(.1996) (BM96). 

To be able to discuss the second possibility laid out 
above, we use the same dynamical collapse model to con- 
struct barriers corresponding to filaments and sheet forma- 
tio n. Such barriers w ere first presented in a recent paper 
bv lShen et alJ <|2006F) . who used them, among other things, 
to provide interesting analytic estimates of the mass frac- 
tions in pancakes, filaments and haloes at any given time. 
Although they used the average shape approximation, some- 
thing we will avoid here, their barriers could in principle also 
have been used for parts of the discussion in this paper. An 
alternative approach using the Zel'dovich approximatio n to 
discus s pancake formation was recently suggested by iLed 
J2006D . 

Our paper is laid out as follows. We first review the 
standard implementation of the excursion set method with 
flat and moving barriers and demonstrate why this is un- 
able to reproduce the correlations discussed above. We then 
explain the generalisation of CL01 to a random walk in the 
deformation tensor. We show how choosing a particular lin- 
ear combination of the eigenvalues allows the behaviour of 
the ellipsoidal collapse model to be successfully captured by 
barrier functions for haloes, filaments and sheets. We then 
run large samples of Monte Carlo realisations and analyse 
the results. We will conclude that it is possible to under- 
stand the clustering dependence on formation history if we 
make the plausible assumption that the latter depends on 
the masses of progenitor pancakes and filaments. 



2 THE EXCURSION SET APPROACH 

2.1 A random walk with a flat barrier 

The spherical collapse model identifies a linear theory den- 
sity contrast, S c (z) = 5 C /D(z), extrapolated to present 
epoch, whic h corresponds to collaps e at redshift z. Here 
S c w 1.686 JPress fc Schechterl 11974) and D(z) is the lin- 
ear theory growth function normalised to 1 at z = 0. This 
equals 1/(1 -j- z) in the special case of an Einstein de Sitter 
universe. 

The traditional excursion set approach entails picking a 
random position (particle) and smoothing the linear density 
field, S, over ever smaller spheres around this point until this 
spherical collapse criterion is fulfilled (8 ^ S c (z)). 

The density perturbation smoothed over a scale R is 
given by 

8r(x) = J 5(x')W R (x - x')d 3 x', (1) 

where Wr(x) is the scale-dependent smoothing filter. The 
variance of the density field on scale R is given by 

a 2 (R) = (\5 R (x)\ 2 ) = JL J p(k)W 2 R (k)k 2 dk, (2) 

where W R (k) is the fourier transform of the filter function. 
We are in principle free to choose this filter function, and the 
top-hat filter is frequently used for visualisation purposes, 
although computationally it is awkward since the added 
modes are correlated as we go to smaller scales. If instead a 



Why does the clustering of haloes depend on their formation history? 3 



hard k— space filter is used (with an upper k— space cutoff 
k c — 1/-R), the modes added by going to a smaller scale 
are statistically independent of the modes already included. 
Thus the density contrast performs a Markov random walk, 
and the entire problem can be described by such a random 
walk with an absorbing constant barrier, where a 2 is the 
pseudo-time variable. 

A particle trajectory will thus gradually carry out a ran- 
dom walk as a 2 increases until the barrier is crossed at some 
particular value of this variance. This value corresponds to a 
smoothing scale R, and the mass of the corresponding halo 
is usually assigned using the top-hat formula (independent 
of the filter actually used for the smoothing) ; 



M 
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For the flat barrier, the problem of a Markov random 
walk with an absorbing barrier is analyticall y tractable. Th e 
first crossing distribution is given by jChandrasekarlll943t) 



^ /M =fe) 1/2exp H 



(4) 



where v = S 2 /a 2 . 

The mass functions n(m, z) of dark haloes are obtained 
from the density distribution function of up-crossings f{v), 
through 



vf{i>) = rn 



2 n(m, z) din m 
p d\nu 



(5) 



It is important to realise that only the first crossing of 
a barrier is of any importance. Any subsequent crossings of 
the same barrier are irrelevant since collapse has already 
occurred at this redshift on a larger scale. However, an ear- 
lier redshift corresponds to a higher barrier, and the first 
up-crossing of this higher barrier corresponds to the largest 
scale on which collapse has occurred at the earlier redshift. 

The conditional mass functions, defined as the distribu- 
tions of progenitor mass at an earlier redshift for given final 
halo mass, are related to this two barrier problem. In the 
flat barrier case they are given by a simple extension of the 
unconditional mass function. 

Although the excursion set approach cannot tell us 
about the substructure of a halo at any particular redshift, it 
does tell us something about how mass was added to the halo 
as it grew. This is possible since we can follow the mass of 
the halo's main pro genitor with redshift , and thereby study 
its merger history llLacev fc Colelll993D . This allows us to 
identify the formation redshift Zf as the redshift at which 
we can first identify a progenitor with more than half the 
final mass. 



2.2 The moving barrier 

ISheth et all (1200 ill (SMTPl! were the first to use the ellip- 
soidal collapse model of BM96 in the excursion set approach. 
In this model collapse depends not just on the density field, 
but on eigenvalues of the tidal shear tensor. However, rather 
than drawing realisations of the shear tensor, and thereby 
including full shape information, SMT01 used expectation 
values of the shape parameters as a function of a. They 
were then able to construct a collapse barrier whose value 
depended on the variance, a 2 , thus on mass. Since variance 
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Figure 1. Example of a 1-dimensional random walk (jagged line) 
together with the average barriers for sheets, filaments and haloes 
at Z = 0. The variance at which first up-crossing occurs deter- 
mines the mass of each structure element in which this specific 
particle is embedded. It is clear that each halo identified at z = 
is contained in a larger mass filament which is in turn contained in 
an even larger pancake. The environment density of the halo de- 
pends on values of 8 on yet larger scales, nearer to a = 0, whereas, 
in the usual formulation, its formation history is determined by 
the walk's structure on small scales to the right of halo crossing. 



is the time-like variable, this has been labelled a "moving" 
barrier. 

With this shortcut, the approach, with a sharp fc-space 
filter, was highly successful in reproducing the shapes of 
the halo mass functions seen in N-body simulations. SMT01 
provided a convenient analytic approximation to the first- 
crossing distribution corresponding to their moving barrier, 
but the symmetry which allows for compact representation 
of the conditional mass functions in the flat barrier case is 
lost. 

As an extension to this work lShen et alj fcOOtJ) recently 
provided moving barriers which can represent collapse along 
one, two and three axes, corresponding to the formation of 
sheets, filaments and haloes, an effort which we replicate and 
extend later in this paper. 



2.3 Correlations between clustering and formation 

It is easy to see why the usual excursion set approach pre- 
cludes correlations between large-scale e nvironment a nd for- 
mation history for haloes of given mass iWhitelll996f) . Take 
the walk in figure as an example. The walk crosses the 
z — halo barrier at a 2 = 4.15, and is so is identified as 
belonging to a halo of mass M(<7h). The density of the large- 
scale environment for this particle can be identified using the 
values of 5 near the origin, a 2 <C a\ (and S <C S c ), whereas 
its formation history (at least, according to the usual def- 
initions) depends on the behaviour of the random walk to 
the right of the barrier crossing at a\ = 4.15. Since, by 
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rier increases with a 2 , the filament barrier is constant, and 
the pancake barrier actually decreases as we go to smaller 
scales. Since going to higher redshift simply scales the barri- 
ers up and to the right, it is possible, even probable for small 
masses, that the mass of a progenitor pancake at some early 
time is larger than the final mass of the halo. Figure shows 
an example of such an occurrence, where the progenitor of a 
z — halo was part of a larger mass pancake at z = 1. The 
idea that halo formation history could somehow depend on 
a different barrier than halo mass, is crucial, since there is 
every reason to expect a clear correlation between properties 
defined by the behaviour of the random walk on the same 
side of the barrier crossing which defines halo mass. 

We have identified two independent extensions which 
can, in principle, accommodate the kind of correlation be- 
tween environment and formation history that we seek. To 
examine whether these extensions can provide a quantitative 
explanation of the numerical results, we need to implement 
them in a generalised excursion set approach. This is the 
subject of the next section. 



Figure 2. The random walk in figure with the z = halo 
barrier and 2 = 1 sheet barrier. This trajectory corresponds to 
a halo whose z = 1 progenitor was part of a pancake with mass 
exceeding the z = halo mass. Assuming the assembly history of 
the final halo to be influenced by the mass of this early pancake 
provides a possible explanation for the dependence of clustering 
on formation history, since both the large-scale density and the 
mass of the early pancake correspond to points on the random 
walk to the left of the barrier crossing which defines the halo's 
z = mass. 



assumption, all haloes of given mass cross the z = halo 
collapse barrier at the same variance, the Markov nature of 
the random walk prevents any correlation between the two 
sides, between environment and formation history. In other 
words, old haloes should cluster in exactly the same fashion 
as young ones. 

In the following we discuss two ways in which the excur- 
sion set approach can be extended to accommodate correla- 
tions between clustering and formation history. In the first, 
a multi-dimensional implementation of the random walk is 
co mbined with a non-sp herical collapse barrier, as pioneered 
bv lChiueh fc Led (1200 ill . This implies that halo collapse at 
given mass is associated with a single constraint on a set 
of variables, in their case {8, r}. Formulated differently, the 
effective barrier for S depends on the value of the shape pa- 
rameter r. At barrier crossing 5 and r can thus have different 
values for different halos (though they are perfectly corre- 
lated) and these different values can, in principle, correspond 
to different distributions of both formation history and large- 
scale environment. A correlation between environment and 
formation history thus becomes possible. 

A second possibility is to stipulate that formation his- 
tory of a halo somehow depends not only on its progenitor 
haloes, but on progenitor pancakes and filaments as well. 
These structures are characterized by collapse along one or 
two axes, respectively, and can be represented in the ex- 
cursion set a pproach by differe nt collapse barriers, as first 
suggested bv lShen et alJ (120061) . Figure Q compares the av- 
erage shapes of such barriers at z — 0. The general trend is 
clearly very different as we go to smaller scales; the halo bar- 



3 THE GENERALIZED EXCURSION SET 
METHOD 

3.1 The 6D Random Walk 

The general non-spherical excursion set approach with a 
hard fc-space filter involves a random walk in the deforma- 
tion tensor (the 2nd derivative tensor of the peculiar gravi- 
tational potential), 



(6) 



rather than simply in the density contrast alone (which is 
the trace of this tensor). CL01 devised an algorithm for this: 
The deformation tensor d can be simulated by drawing 
six independent Gaussian variables {yi-.ye} with dispersion 
ctq, and using the following linear transformation. 



d 2 2 = 



yi + 



3 l 

i 5 y2 + 7E y:i 



- 3 \yi 
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15 



di 
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(7) 

(8) 

(9) 
(10) 

(11) 
(12) 



This transformation satisf ies the correlat i ons o f the defor- 
mation tensor as shown in lBardeen et all l)l98f^ . 

The random walk proceeds by drawing new values from 
the Gaussian distribution, and adding to yi, until we have a 
walk of N steps. To each step n can be assigned a dispersion 

2 2 

a = na Q . 

For each step of the walk the deformation tensor can 
be diagonalised to find the three eigenvalues {Ai,A2,A3}. 
This set of eigenvalues, along with a non-spherical collapse 
model, determine whether collapse has occurred by a given 
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redshift. Although the eigenvalues could be drawn directly 
from t he probability distribution function of lDoroshkevich| 

15 3 1 / 3/ 2 15/2 \ 

p(Ai, A2, A3) = —-^exp^^ + ^j (13) 

x(Ax-Aa)(Aa-Aa)(Ai-Aa), (14) 

the above approach is significantly faster. 

CL01 chose to avoid the diagonalisation process alto- 
gether by expressing the collapse criterion in terms of S and 
of a variable r 2 directly calculable from the y^s through two 
rotational invariants. This speeded up their calculations and 
allowed a nice visualisation of the problem. We discovered 
that this simplification results in a suboptimal representa- 
tion of the fully 3-dimensional ellipsoidal collapse model, and 
we choose instead to use diagonalisation of the deformation 
tensor to obtain the the eigenvalues Ai,A2,A3. The change 
of variables Ai, A2, A3 — > 5,v,w given by 

5 = A1 + A2 + A3 (15) 
v = -Ax + Aa (16) 
w = -Ai-A 2 + 2A 3 (17) 

allows simple and accurate modelling of the behaviour of the 
ellipsoidal collapse model. This choice is discussed in detail 
in Appendix IbI 

3.2 Collapse Barriers 

The spherical collapse barrier is found by applying the spher- 
ical collapse model and identifying the value of the initial 
density contrast, linearly extrapolated to present time, which 
corresponds to collapse at redshift z. It thus represents a 
constant (flat) barrier for the random walk. In the ellip- 
soidal collapse model, collapse is determined by the shape 
of the initial patch as well as its overdensity and this results 
in a barrier which depends on all 3 eigenvalues Ai rather 
than just on their trace 5. 

To fully capture this dependency we use the following 
barrier shape to represent halo collapse. 

S B (v, w, z) = -A_ [1 + Ql (vD(z)) a2 + q 3 [wD{z)) a *] (18) 

where S c = 1.686 is the critical value for spherical collapse 
at the present epoch, D(z) is the linear theory growth factor 
normalised to 1 at z = 0, and 

Qi « 0.2809, q 2 w 1.3557, a 3 « 0.070, a 4 « 1.41205 (19) 

The coefficients here were found by fitting to the col- 
lapse values of 8000 walks (stepsize Act 2 = 0.025) for which 
we used the full numerical solution to the ellipsoidal collapse 
model of BM96 (see appendix). This interpolation formula 
provides a near perfect fit to the true barrier obtained using 
the numerical solution, and it allows for a thorough investi- 
gation of the limitations and accuracy of the excursion set 
approach. We want to stress that our barrier was not found 
by fitting to N-body mass functions, but rather by requiring 
a n accurate repre s entat ion of the ellipsoidal collapse model 
of iBond fe Mverd lll996Tl . 

We also present barrier functions representing collapse 
along one and two axes. For sheets we have 



Ss(v, w, z) = [1 - ai (vD(z)r - MwD(z)) a *] (20) 

D{z) 

where on = 0.3748, « 2 = 0.2399, « 3 = 0.003237, « 4 = 
2.4187. The formula for the filaments is slightly more com- 
plicated since it needs to reproduce halo-like behaviour for 
two large axes, and sheet-like behaviour for two short axes. 

Sf(v, w, z) = [1 + a lV D(z) (1 + a 2 (xD(z)) a3 )] (21) 

where ai = 0.1104, a 2 = 0.01641, a 3 = 1.463, and x = 
w + 3v,y — w — 3v. This change of variables is meaning- 
ful, since the condition of perfect triaxiality, for which the 
behaviou r of the second ax is is nearly identical to the spher- 
ical case llShen et alj|2006l) . translates to w = 3i>. Using the 
expectation values found in Appendix [B] we c an translate 
these barriers into "moving barriers" similar to IShen et alJ 
(2006). These average barriers are the ones plotted in figure 
□ 



4 MONTE CARLO SIMULATIONS 

To address the questions posed in section |2] we simulate a 
large ensemble (~ 10 6 ) of random walks. For each simulated 
trajectory we output the density contrast on large scales 
along with the first up-crossings (actually the correspond- 
ing masses) of the halo, filament and sheet barriers fsec. l3.2t 
for a range of different redshifts. We also store the nominal 
formation redshift for the z — halo, as well as the values 
of the auxiliary variables v and w at collapse. In the follow- 
ing we analyse the results and we try to understand their 
implications. 

4.1 First crossing distributions 

Figure @ shows the simulated first crossing distributions 
for the halo barrier introduced in section 13.21 (squares) 
and for the spherical collapse barrier (triangles). These are 
compared with the corres ponding analytic predictions from 
iPress fc SchechteJ lll974fi and SMT01. It should come as 
no surprise that the distribution for our ellipsoidal collapse 
model is well fit by the SMT01 analytic formula, since this 
was obtained by fitting to a "moving barrier" version of 
our approach. Note that their final mass function contained 
an additional scaling parameter a, which allowed for better 
agreement with simulations. This degree of freedom is ab- 
sent in our approach, though one might argue that it should, 
perhaps, be reintroduced when comparing our predictions to 
numerical results obtained with any particular (and to some 
extent arbitrary) definition of the boundary and thus the 
mass of a halo. 

4.2 Correlations between large and small scales 

As laid out in section l2~3l a multi-dimensional random walk 
with a shape-dependent barrier allows, in principle, for cor- 
relations between the environment density on large scale and 
the point at which the collapse barrier is crossed. 

In the following we use our Monte Carlo sample of ran- 
dom walks to quantify this connection between large-scale 
overdensity, point of barrier crossing and, potentially, halo 
formation redshift. 
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Figure 3. First crossing distributions obtained using the spheri- 
cal collapse barrier (red squares) and the ellipsoidal collapse bar- 
rier presented in this paper (black triangles). The solid line is 
the analytical Press-Schechter solution, and the dotted line is the 
Sheth-Mo-Tormen approximation. 

4-2.1 Environment density and point of barrier crossing 

From our ensemble of walks we pick out those in the up- 
per and lower 10% tails of the environment overdensity dis- 
tribution (estimated as S(a 2 — 0.5)). This represents, in a 
statistical sense, our most and least clustered haloes. We 
then compare the distributions of v/a and w/a at barrier 
crossing for the two sets. Concentrating on ratios removes 
the obvious dependence on a, providing shape parameters 
which allow us to compare positions of barrier crossing in a 
way which does not depend on overall barrier scale. 

For the ellipsoidal collapse barrier we find a clear albeit 
small correlation between environment and point of crossing. 
Both w and v are on average higher for walks with high 
environment density. The average value of v/a is roughly 
7 — 8% larger for the walks with high environment density 
than for walks in the low density tail. The dependence on 
w is somewhat less pronounced (5 — 6% difference between 
the two tails). As expected no such effects are found when 
we use the spherical collapse boundary. 

It is worth noting that this effect is stronger for low 
values of a 2 at barrier crossing (i.e. for massive haloes) 
than for less massive objects. This can be understood by 
recalling that the random walk is a diffusion process, so 
that correlations with values at "early times" (i.e. large 
scales, small variance values) are gradually washed out 
as "time/distance" increases (i.e. as variance increases). 
Throughout this work we define the environment by the 
overdensity at a fixed scale a 2 = 0.5, roughly corresponding 
to R ~ 10h~ Mpc in the concordance ACDM cosmology we 
are adopting. With this definition it is not unexpected that 
environment correlates less with properties at barrier cross- 
ing as we go to smaller scales (larger a 2 and smaller masses) . 
This gives a first hint that the 6D excursion se t approach 
will be unable to explain the effects measured bv lGao et alJ 
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Figure 4. Histogram of e values at first crossing for the 10% 
highest (black) and 10% lowest (red) tails of environment density 
(estimated at a 2 = 0.5) for the ellipsoidal collapse barrier. This 
plot refers to haloes more massive than the characteristic mass 
M*. Although the effect is slight, the difference between the dis- 
tributions is statistically highly significant. The values of e are on 
average ~ 10% larger at the point of first crossing for walks in the 
high environment density tail than for walks in the corresponding 
low density tail. 

(2005), since these are largest for low mass haloes, and are 
virtually non-existent for the most massive haloes. 



4-2.2 Clustering and Patch Shapes 

It is perhaps more interesting to study environmental effects 
on collapse properties we are familiar with. For the purpose 
of this analysis we estimate the ellipticity and prolateness of 
the initial density field associated with a halo by 

e = (-Ai + A 3 )/<t (22) 

and 

p = (Ai - 2A 2 + A 3 )/cr (23) 

We find that there is a weak correlation between en- 
vironment density and the pre-collapse ellipticity of haloes 
(see eqn. 1221 for the relevant definitions). Again we see that 
the effect fades as we go to smaller masses. Figure (0J shows 
histograms of e values for the 10% tails of the distribution in 
environment density for haloes with M > M* . A correspond- 
ing plot for haloes with M < 10 _2 M» shows no discernible 
difference between the two distributions. Roughly speaking, 
for massive haloes, those with the 10% highest (lowest) envi- 
ronment densities have ellipticities which are systematically 
5% higher (lower) than those of typical haloes. Systematic 
effects with prolateness are also clear, albeit smaller, indi- 
cating that haloes in dense environments have a wider than 
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typical range of prolateness values - in agreement with the 
results for the ellipticities. 

These results show shows that haloes in dense regions 
originate from more elliptical initial patches than less clus- 
tered objects of similar mass. It is tempting to extrapolate 
this conclusion to apply to the actual ellipticity of the col- 
lapsed, quasi-equilibrium halo. This is speculative, however, 
since we do not know how (or even if) the final shape of a 
halo is related to the shape of its initial patch. 

4-2.3 Environment Density and Formation Redshift 

In sections 14.2.11 and I4.2.2H we established that there is a 
small but clear correlation between environment overden- 
sity and the point of barrier crossing, thus with the shape 
of the collapsing patch. The point of barrier crossing clearly 
influences the properties of the random walk at larger vari- 
ance, hence on smaller mass scales, so we may expect a cor- 
relation between environment density and formation red- 
shift. T his is precis e ly wh at was detected in N-body simula- 
tions bv lGao et alJ (120051) . They show that low-mass haloes 
with high formation redshifts are significantly more clus- 
tered than similar mass but late-forming objects. Hence they 
sit in denser environments. 

We want to compare our results as directly as possible 
with the results of GSW05. They look at a sample of haloes 
with particle numbers in the range 100-200. Each particle 
has a mass of 8.6 x lO 8 ft _1 M0, so this corresponds roughly 
to ~ 1 — 2 x 10 11 hT 1 Mq. In terms of the characteristic mass 
M* = 6.15 x 10 12 M Q this is roughly M = 0.14 - 0.28 x M». 

We try to replicate this by taking similar samples from 
our random walks. We use an effective power spectrum in- 
dex of n — —2, and look at two mass bins, with masses 
in the ranges M = [M*,2M„], and M = [0.02M,, 0.1M,]. 
To improve statistics we use larger bins than GSW05 and 
we compare environment overdensity for haloes in the 20% 
(rather than 10%) tails of the distribution of formation red- 
shift. 

Despite these attempts at improving our statistical 
leverage, we are unable to find any effect whatsoever. We 
strongly exclude anything of the magnitude seen by GSW05. 
Our results therefore appear to rule out this idea as a po- 
tential model for the effects seen in the simulations. 

This should not take us entirely by surprise, since the 
results of the previous section lead us to expect that any 
effect should be larger for the high-mass than for low-mass 
haloes. This is opposite to the numerical result. We conclude 
that although a formal dependence should be present, it is 
far too small to measure and cannot explain the results of 
GSW05 and CGW06. 



4.3 Clustering and progenitor pancakes and 
filaments 

Since our multi-dimensional extension of the excursion set 
approach is unable to reproduce the required correlation be- 
tween clustering and formation history, we now turn to the 
other possibility laid out in section 12.31 namely that the 
correlation may reflect a dependence of formation history 
on the properties of progenitor pancakes and filaments. In 
particular we study the correlation between the masses of 




M / M" 



Figure 5. Average density at large scale (i?o ~ 10h -1 Mpc) 
around a z = halo, as a function of halo mass (black curve) 
together with the same quantity for halo subsets defined so that 
the progenitor pancakes at z = 2 lie in the upper (blue curve) 
and lower (red curve) 10% tails of the pancake mass distribution. 



these progenitor structures and the large-scale environment 
of the z = halo. 

Figure |K| plots average large-scale (~ 10/i _1 Mpc) envi- 
ronment overdensity as a function of halo mass, and com- 
pares this with the values found using only those haloes 
whose progenitor pancakes at z = 2 lie in the upper or in 
the lower 10% tail of the pancake mass distribution. There is 
a strong correlation between this large-scale overdensity and 
the mass of the progenitor pancake, with haloes in denser 
environments tending to have more massive progenitor pan- 
cakes. The strength of this effect is greatest for low-mass 
haloes. 

Fig. HJ shows a similar plot of large-scale environment 
overdensity as a function of halo mass, but now sorting 
haloes by progenitor filament mass at z = 2. A dependence 
is again clear, but is much weaker than was the case for pro- 
genitor pancakes. The smallness of the effect in this case is 
due to the fact that very few haloes were embedded in higher 
mass filaments at z = 2. Only for these haloes is any cor- 
relation between environment and progenitor filament mass 
to be expected, as explained in sec. 12.31 . 

We have established a clear correlation between envi- 
ronment overdensity at large scales and the mass of progen- 
itor pancakes and filaments. If the formation history of a 
halo is somehow affected by the size of its progenitor pan- 
cakes and filaments, this correlation will induce a correlation 
between clustering and formation history. Intriguingly the 
effect is largest for low-mass haloes and becomes quite weak 
for the most massive haloes. This is precisely the behaviour 
found by GSW05 and CGW06, where the measured effects 
are only strong for haloes less massive than ~ 1O 13 M0. In 
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Figure 6. Average density at large scale (i?o ~ 10h~ 1 Mpc) 
around a z = halo, as a function of halo mass (black curve) 
together with the same quantity for halo subsets defined so that 
the progenitor filaments at z = 2 lie in the upper (blue curve) 
and lower (red curve) 10% tails of the filament mass distribution. 

the schematic treatment of this section, the mass at which 
the dependence becomes small is dependent on whether we 
consider filaments or pancakes, and on the redshift at which 
we consider them to be relevant. 

It is beyond the scope of the current discussion to 
present a detailed model for the dependence of halo for- 
mation history on the properties of progenitor pancakes and 
filaments. Such a connection is intuitively appealing, how- 
ever, and we expect this to be an area where further analytic 
modelling will be very fruitful. 

4.4 Pre-heating by pre-virialization? 

Our barriers for collapse alon g one or two axe s also allow us 
to address a recent claim bv lMo et al.l i2005tl that the ma- 
terial which now resides in low-mass haloes (M < 1O 12 M ) 
was pre-heated by the collapse of larger-scale pancakes (typ- 
ical mass M ~ 5 x lO 12 M fo _1 ) at z ~ 2, and that as a result 
the gas component of many of these haloes failed to cool and 
form stars. The analytic argument underlying this claim was 
based on a simplified calculation which specifically failed to 
include the conditional distribution of the masses of progen- 
itor pancakes for haloes of given z = mass. We are now 
well equipped to carry out a more precise calculation of the 
relevant quantities. 

To examine this claim we pick out haloes in narrow 
mass ranges around M = 1O 12 M , M = 10 11 M and 
M = 10 10 M o , and we look at the mass distribution of their 
progenitor pancakes and filaments at redshifts 1, 2, and 3. 
These cumulative distributions are shown in figs. J7J (pan- 
cakes) and (JHJ (filaments). They show that many low-mass 
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Figure 7. Cumulative distribution of progenitor pancake masses 
at several redshifts for haloes of z = mass 10 12 Mq (solid), 
10 n M Q (dashed) and 1O 1O M (dotted). At z = 2 only ~ 10% of 
10 12 Mp, haloes had progenitor pancakes of the size suggested by 
iMo et alJ i2005ly For 10 n Af Q and 1O 1O M haloes this fraction 
is significantly smaller. Note that almost all haloes in these mass 
ranges were embedded in larger mass pancakes at z = 1. At z = 3 
this fraction is 30% for the 1O 12 M0 haloes but is still large (~ 
80%) for the smaller mass ranges. 
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Figure 8. Cumulative distribution of progenitor filament masses 
at several redshifts for haloes of z = mass 1O 12 M0 (solid), 
1O U M (dashed) and 1O 1O M (dotted). 
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haloes were indeed embedded in larger mass str uctures at 
redshi fts z ~ 1 — 2, in qualitative agreement with lMo et alJ 

Quantitatively however, the m echanism looks less 
promising. Central to the argument of iMo et alJ 1120051) was 
the assumption that most low-mass haloes were embedded 
in pancakes of mass M ~ 5 x 10 12 Af Q at z ~ 2. This was 
required to ensure sufficient preheating of gas to affect later 
condensation and star formation. From the middle graph in 
fig. J7J it is clear that typical IO^Mq haloes had progeni- 
tor pancakes of mass M ~ 2 x 10 Mq at Z = 2, and only a 
fe w percent have pr ogenitor pancakes as massive as assumed 
bv lMo et alJ (120051) . The situation is even worse for typical 
10 10 M Q haloes which were embedded in pancakes of mass 
M ~ 7 x 10 10 Mq at z = 2, nearly two orders of magnitude 
less than the required value. 

Our results contradict those of IMo et alJ 1120051) in two 
ways. According to our analysis, progenitor pancake mass is 
not independent of halo mass in the way they suggest, and 
even for relatively massive "dwarf galaxy" haloes, the pro- 
genitor pancakes are typically too small to generate the tem- 
peratures required for suppression of later cooling. Taking 
bot h these discrepa ncies into account, it seems unlikely that 
the lMo et alJ (120051) mechanism can account for the shallow 
faint end slope of the galaxy luminosity function and the rel- 
atively small number of luminous satellites around galaxies 
like the Milky Way. 



5 CONCLUSIONS 

We have used the excursion set formalism to discuss the cor- 
relation between environment and formation history found 
for haloes in large-scale simulations of cosmic structure for- 
mation. We implemented a multi-dimensional generalisation 
of the excursion set approach which assu mes halo forma- 
tion t o follow the triaxial collapse model of iBond fc Mversl 
(1996) and which does not require further simpl i fying as- 
sumptions of the kind introduced bv lSheth et alJ (1200 if) to 
reduce the system to a 1-dimensional barrier crossing prob- 
lem. A correlation between halo formatio n history and halo 
environment is expected in this model (IWang et alJ 1200111) 
but we have shown here that it is very weak. We did, how- 
ever, demonstrate that there is a measurable correlation be- 
tween environment density and the pre-collapse shape of a 
halo. 

It appears that haloes in dense regions have slightly 
more extreme axial ratios than similar mass objects in un- 
derdense environments. This could potentially be studied 
with N-body simulations, although we caution that our def- 
inition of ellipticity, eqn. ]22L while mathematically conve- 
nient, is not simply related to the measurable shape of non- 
linear, quasi-equilibrium haloes. To address this issue prop- 
erly, we suggest using la rge N-body simulatio ns such as the 
Millennium Simulation dSpringel et aljfeoOSl) . It would be 
feasible to follow particles which are part of z — haloes 
back to the initial density field. With a large enough sam- 
ple of haloes one could then determine how the final shapes 
of halos are related to the properties of their initial, linear 
"patches" . 

With straightforward definitions of environment den- 
sity and halo formation time our multi-dimensional random 



walks are inconsistent with the environment-formation his- 
tory correlations found by GSW05 and CGW06. We have, 
however, discovered a different mechanism which can induce 
a dependence of some aspects of halo formation on halo en- 
vironment. Crucial to our analysis of this mechanism are 
the collapse barriers for haloes, filaments and sheets intro- 
duced in section 13.21 . We showed that many low-mass z = 
haloes were embedded in larger-mass progenitor pancakes at 
redshifts z ~ 1 — 2, and that the mass of these progenitors 
correlates with the density of the environment surrounding 
the z — haloes. Thus, if the formation history of z — 
haloes is affected by the mass of their high-redshift progen- 
itor pancakes, we can expect an induced dependence of for- 
mation history on environment. Although we do not attempt 
to model this, we show that sorting haloes by the mass of 
their progenitor pancakes at z ~ 1 to 3 gives a variation 
of mean environment with pancake mass which is of similar 
size and has a similar dependence on halo mass to the vari- 
ation found in simulations when haloes are ranked by their 
formation time. Future efforts should result in a more de- 
tailed model for the connection between formation history 
and the masses of progenit or pancakes and f ilaments. 

A very recent paper bv lWang et alJ i2006f) suggests that 
the age dependence of clustering observed by GSW05 is re- 
ally caused by mass stripping due to the large scale tidal 
field, which makes small old haloes appear less massive than 
their initial density field would suggest. This might provide 
an intuitive explanation for the effect, and is an example of 
an effect which could be connected to our scenario. 

By studying conditional mass functions for the sheets 
and filaments in which the material of low-mass z — haloes 
was embedded at early times, we found typical masses for 
progenitor pancak es at z = 2 which are significantly smaller 
than suggested bv lMo et al.l (120051) . In addition, these pro- 
genitor masses depend strongly on the mass of the final halo. 
These results suggest that preheating due to the collapse of 
progenitor pancakes is insufficient to suppress later cooling 
and condensation of gas. It is nevertheless interesting that 
the sheet mass function is more strongly peaked than its 
halo counterpart, so that the majority of low-mass haloes 
were indeed embedded in more massive pancakes at z ~ 2. 

In summary we argue that this work may provide an 
interesting and intuitive route to understanding how envi- 
ronment influences halo formation. It is important that an- 
alytic methods are refined to account at least qualitatively 
for the relatively strong effects seen in N-body simulations, 
and our results provide a first step in this direction. To build 
on the approach presented here it will necessary to model in 
more detail how the structural elements of the cosmic web 
influence the formation of haloes within them. We hope our 
results will motivate such efforts. 
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APPENDIX A: 
MODEL 



THE ELLIPSOIDAL COLLAPSE 



For convenience we here review the most important features 
of the ellipsoidal collapse model of BM96. 

Starting from a homogeneous ellipsoidal patch with 
density contrast 5, the general equation that governs the 
collapse of an initially spherical surface along principal axis 
k is 



dt 2 



-AirGpA k 



1+5 , K 



(Al) 



where b' k 5/2 and A^ denote the interior and exterior tidal 
forces respectively. The interior tidal forces can be calculated 
for a constant density patch from the potential theory of 
homogeneous ellipsoids (cf. Binney & Tremaine 1987), and 
are given by: 



dr 



(4 + r) 



OJ+T 



The linear approximation for external tides is: 



A extk — 



(A2) 



(A3) 



where the A^s are the initial eigenvalues of the deformation 
tensor. The evolving quantity Ak relates to the boundary 
radius of the Lagrangian sphere, R through Ak(t) = a,k(t)R, 
so that the a k can be thought of as the anisotropic general- 
isation of the standard cosmic expansion factor. We should 
note here that these equations are valid only in cosmological 
models with vanishing cosmological constant - a discussion 
of the modifications for a non-zero cosmological constant 
can be found in BM96. 



The initial conditions are set by the Zel'dovich approx- 
imation 



A k (ti) = di [1 - \ k (ti)] 

A k {U) = H(ti)A k (ti)-aiH(ti)X k (ti 



(A4) 
(A5) 



The differential equations Al are evolved for all axes k 
until a value of Ak = 0.177a has been reached, after which 
the individual Ak are held constant. The factor of 0.177 re- 
sults from the requirement that the virial density contrast 
of 179 obtained from spherical top hat calculations is repro- 
duced. 



APPENDIX B: INVARIANTS AND VARIABLES 

The three rotational invariants of the 3x3 deformation ten- 
sor are 



11 = Ai + A2 + A3 

12 = A1A2 + A2A3 + A3A1 

13 = A1A2A3. 



(Bl) 
(B2) 
(B3) 



They can all be easily expressed in terms of the yiS. It would 
therefore be useful if we could express our barrier variables 
in terms of these invariants. 

The first barrier variable is naturally the density con- 
trast, 8 = Ii =j/i. The second variable chosen by CL01, is 
a combination of the first two invariants 



= 3/2/ 1 2 -2/ 2 = i^(A i -A J ) 2 



(B4) 



which satisfies the condition that it is independent of y\. 
This is important both since we need to fit the density con- 
trast at collapse to the barrier parameter, and because it 
means that it vanishes in the spherical case. (The spherical 
case is equivalent to yi^i =0.) 

For the third barrier variable a similarly convenient 
combination of the three invariants would be 



= - (2/? 
9 v 1 



97i/ 2 + 27/3 



(B5) 



which is also independent of y\ and thus also vanishes in the 
spherical case. 



The two variables 
from the ys as follows: 



and u may be obtained directly 



(2/2 + 2/3 + yi + yt 



Ve) 



(B6) 



*>ylyi 



.'/.-//:! 



+ 6\/5' 



2/6 2/3 



1/42/3 + 6\/52/22/3 



^(-3V5V3, 
+3VEV3y 4 y2 — 3\/52/4 

+2 - V5y! - 6V5V3J/4J/52/6) (B7) 

_where we have used the explicit formulae in lLin et alJ 
)(LCL02). We thus have two variables which are ex- 
pressed in terms of the rotational invariants and which are 
independent of the density contrast. The question is whether 
a barrier in these two variables is capable of capturing the 
behaviour of the ellipsoidal collapse model, and here we 
stumble upon a problem: 

Although the transformation involved is one-to-one it 
does not guarantee that it will allow one variable to be a 
single valued function of the two others. Indeed it is easy 
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to show that transforming from eigenvalue space to 5, r 2 , u 3 
occasionally yields identical r and u values for different 8 c s 
- thereby preventing an expression of 5 C as a single valued 
function of the other two variables. 

It proves significantly more accurate to express the bar- 
rier variables as suitable orthogonal linear combinations of 
the eigenvalues as follows 

S = Ai + A 2 + A 3 (B8) 
v = -A1 + A2 (B9) 
w = -A!-A 2 + 2A 3 . (BIO) 

The variables v and w both vanish in the spherical case, are 
independent of S and of each other, and are always positive 
for ordered eigenvalues Ai < A2 < A3. Since we prefer to 
work with a barrier which reproduces the ellipsoidal collapse 
model of BM96 as well as possible, we use this set of variables 
throughout our work. 

The joint probability distribution can be found by ma- 
nipulation of eqn. 114H 

, s 1 v 1 15 3 1 / 1 ,. 3 15 2 5 2 

P (S,v,wW) = 3vT^;exP{-w {6 + ~ V + l W 

x(w 2 ~v 2 )v, (Bll) 
where the limits are 
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(B12) 





< v < 
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(B13) 





< w < 


00 


(B14) 



The expectation values of v and w for a given a are 
(«) = -S= (B15) 

V 1U7T 

{w) = -^L=. (Bie) 

V 1U7T 

These are the values used to obtained the average barrier 
functions in figure Q in the main text. 
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